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Abstract. Recently much effort has been made towards the introduction of non- 
Hermitian random matrix models respecting PT-symmetry. Here we show that 
there is a one-to-one correspondence between complex PT-symmetric matrices 
and split-complex and split-quaternionic versions of Hermitian matrices. We 
introduce two new random matrix ensembles of (a) Gaussian split-complex 
Hermitian, and (b) Gaussian split-quaternionic Hermitian matrices, of arbitrary 
sizes. We conjecture that these ensembles represent universality classes for PT- 
symmetric matrices. For the case of 2 X 2 matrices we derive analytic expressions 
for the joint probability distributions of the eigenvalues, the one-level densities 
and the level spacings in the case of real eigenvalues. 


In recent years there has been a surge of research interest in PT-symmetric 
quantum theories, accompanied by a multitude of experimental applications and 
realisations [lT7 


For finite-dimensional systems represented by matrices, PT- 


symmetry is equivalent to the reality of the characteristic polynomial 18 . That 


is, PT-symmetric matrices have either real or complex conjugate eigenvalues. Their 
eigenvectors are orthogonal with respect to a suitably defined CPT inner product 
[T^ . It has recently been conjectured that PT-symmetry is closely related to 
split-quaternionic extensions of quantum theory 20 21 . Here we show that split- 


quaternionic extensions of Hermitian matrices are indeed a natural representation of 
PT-symmetric matrices. This equivalence allows us to introduce new PT-symmetric 
random matrix ensembles. 

In conventional quantum systems, random matrices play an important role due to 
their ability to describe spectral fluctuations in sufficiently complicated systems 22 


In particular, the famous Bohigas-Giannoni-Schmit conjecture states that the spectral 
fluctuations of quantum systems with chaotic classical counterparts are similar to 
those of certain random matrices 23,24. There are three important universality 


classes for Hermitian quantum systems, depending on the time-reversal properties 
of the system, corresponding to the Gaussian orthogonal, unitary, and symplectic 
ensembles 25 . Non-Hermitian random matrix models, on the other hand, whose 


eigenvalues are in general complex, are widely studied, and have applications ranging 
from dissipative quantum systems and scattering theory to quantum chromodynamics 
and references therein). Several attempts towards defining PT- 


26 


(see, e.g., 

symmetric random matrices and identifying universality classes for PT-symmetric 
systems have been made 


27-32 . Most of them are restricted to 2 x 2 matrices, due 


to the lack of a natural parameterisation of larger PT-symmetric matrices. Here we 
introduce the split-complex and split-quaternionic versions of the Gaussian unitary 
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and symplectic ensembles as candidates for representing new universality classes for 
PT-symmetric quantum systems. The new ensembles are invariant under orthogonal 
and unitary transformations, reflecting the invariance of split-complex and split- 
quaternionic Hermitian matrices, respectively. We show that the split-complex case is 
closely related to the well-known real Ginibre ensemble 26 33 , which in itself is PT- 


symmetric. The formulation in terms of split-complex Hermitian matrices might offer 
a new approach towards solving some outstanding problems regarding the spectral 
features of Ginibre matrices. To the best of our knowledge the split-quaternionic 
Hermitian ensemble is not directly related to any previously studied ensemble. 

Let us start with a brief summary of some important properties of split-complex 
and split-quaternionic numbers and matrices that we will need in the following. Split- 
complex numbers 


34 


are numbers of the form z = x -\- ]y, where x,y G H and j is 
the imaginary unit of the algebra, such that = 1. The conjugate of a split-complex 
number is given hy z = x — iy. The “norm” of a split complex number, defined by 
|zp = zz = — y^, can be positive, null, or negative. A split-complex number can 

be represented as a real 2 x 2 matrix 




( 1 ) 


Similarly, a generic split-quaternion p € IH 5 can be written as p = po + ipi +i\P 2 + 

( 2 ) 


kp 3 , where Pi G R, and i,j,k satisfy the relations 


= - 1 , j = k = ijk = -bl. 


The conjugate of a split-quaternion is defined as p = po — ipi — jp 2 — kp 3 . The 
“norm” of a split-quaternion is defined similar to the split-complex case by |pp = 
PP = Po + p\ ~P 2 ~Psy again can be positive, null, or negative. A split-quaternion 
can be represented as a complex 2 x 2 matrix 


f Po + m P2 + m \ 
\ P2- m Po - ipi j 


( 3 ) 


In analogy to the case of quaternionic vector 
an (indefinite) inner product for two vectors u,v G 
components as 


spaces 36 we can define 
with split-quaternionic 


N 

{u,v) = '^U„Vn. (4) 

n—1 


The adjoint of a split-quaternionic matrix A is then defined as {u,Av) = (A^u, u), 
which in terms of the elements means taking the transpose and split-quaternionic 
conjugate, A^ = = (ofej"). We can generalise the concept of Hermiticity to split- 

quaternionic matrices, by referring to matrices with split-quaternionic elements that 
satisfy = H as split-Hermitian. The set of split-Hermitian matrices is invariant 
under standard unitary transformations. The subset of split-complex Hermitian 
matrices is invariant under orthogonal transformations only. The space oi N x N 
split-Hermitian matrices is {2N'^ — Af)-dimensional, and the subspace of split-complex 
Hermitian matrices is A^^-dimensional. 

We can define eigenvalues and eigenvectors of split-Hermitian matrices using the 
complex 2x2 matrix representation in (|^ for the elements. The resulting characteristic 
polynomial is real, thus the eigenvalues are either real or come in complex conjugate 
pairs. Further, the eigenvalues of the 2N x 2N representation are doubly degenerate. 
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in analogy to Kramer’s degeneracy for quaternionic Hermitian matrices 21 22,24 


The split-quaternionic eigenvectors belonging to distinct eigenvalues are orthogonal 
with respect to the inner product (|^. 

Following 18 we define the class of PT-symmetric matrices as the set of complex 
matrices with real characteristic polynomial. Thus, any split-Hermitian matrix can 
be interpreted as a PT-symmetric Hamiltonian. A simple counting argument shows 
that the spaces of split-Hermitian and complex PT-symmetric matrices are in fact 
isomorphic. An NxN complex matrix has 2N'^ real parameters. The coefficients of the 
characteristic polynomial are linearly independent linear functions of the parameters. 
That is, if we condition the matrix to have a real characteristic polynomial, we find 
one constraint for the parameters for the coefficient of each power in the polynomial, 
yielding N linearly independent constraints. That is, a PT-symmetric matrix can be 
parameterised by 2iV^ — N parameters, just as a split-Hermitian matrix of the same 
size, thus they are isomorphic. In particular, the subspace of split-complex Hermitian 
matrices is isomorphic to the space of real PT-symmetric matrices (which is just the 
space of all real matrices). 

Let us now introduce random matrix ensembles of split-complex and split- 
quaternionic Hermitian matrices, as direct generalisations of the standard Gaussian 
unitary and Gaussian symplectic ensembles. That is, we consider split-quaternionic 
and split-complex matrices whose elements are independently distributed normal 
random variables, subject to the constraint of split-Hermiticity. We define the 
probability density function on the space of split-complex Hermitian matrices H as 


¥{H) dH = ( ^ ) " ( ^ ) 


'dP, 


(5) 


i) 11™=! 


Thus, this ensemble is 


where dH = I\m<n dRe(/l mn ) dlm(/i 

n 

invariant under orthogonal transformations of H. 

In the split-quaternionic case, we define the probability distribution on the space 
of matrices H as 

' 2 \.^ / 2 \ \ \ ^ \ — ± I _( 6 ) 


P(P) dH = (-) " ( 2 ^2W(Ar l)^-Tr{HH^+H^H) 


where dH = J([ 


dh^ d/i' dW d/i*" 

m<n^ ^rnn ^'^mn ^'^mn ^'‘'m 


. 1 Vm—I 


dh„ 


and here H^ performs the 


transpose of H and complex conjugation with respect to the imaginary unit i only. 
Hence this ensemble is invariant under unitary transformations of P, reflecting the 
invariance class of split-quaternionic Hermitian matrices. Note that the matrix 
elements of P are by definition statistically independent for both ensembles © and 
Q. That is, these new ensembles fulfil the two defining properties of the standard 
Gaussian ensembles, of independence of the matrix elements and invariance under the 
invariance transformations of the class of matrices on which the ensemble is defined. 
We conjecture that the two split-Hermitian ensembles constitute universality classes 
for PT-symmetric quantum systems. 

We shall now derive analytic expressions for the joint probability density of the 
eigenvalues and the one-level densities for the case of P = 2, following a method 
similar to the one used in 37 for the real Ginibre ensemble. Let us start with the 


split-complex case. Writing 

Ai 

(5-f jq 


P = 


<5-j7 

^2 


(7) 
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with 7 G IR,,j^ = 1 , we have 


^ ^g-(A?+A^+25"+272 


( 8 ) 


The eigenvalues of H are given by 


Ai + A 2 

A 1.2 - ± 




(9) 


'Ai - A2\2 
2 

Clearly these eigenvalues can be real or complex conjugate, depending on the values 
of S and 7 . 

To deduce the joint probability density of the eigenvalues Ai _2 we have to 
transform from A 1 . 2 , 7 , and S to Ai_ 2 ) 7 , and S, and then integrate over 7 and S. 
We have 


Ai,, = ^i±^± 




and thus we find 


P(Ai,2,(5,7) = 


- 52 + 72 , 

1 |Ai - A2|e“('^i+'^2+47^) 


( 10 ) 


( 11 ) 


'(AifAa) -S2+J 

When integrating over 5 and 7 we have to consider that Ai,A 2 , 5,7 are real by 
definition. That is, we integrate over the regions 5^ < ^ + 7 ^, and 7 ^ > 


to find 


P(Ai,A 2 ) = 


e ('^i+'^2)|Ai — A2|erfc(2|Im(Ai_2)|) 

2^/7^ 


( 12 ) 


From the joint probability distribution we can deduce the one-level density of 
the eigenvalues. For this purpose we have to distinguish between the case of real or 
complex conjugate eigenvalues. If Ai,A 2 are real then erfc(O) = 1. Integrating over 
one of the real eigenvalues then yields the one-level density for real eigenvalues 


i?f(A)=i?i(A|AeR) = 


Ae 




-erf(A) 


^- 2 A^ 


(13) 


2 ' 2^ ■ 

We see that there is a non-zero probability to obtain two real eigenvalues, which is 
given by the normalisation of i?®(A), as 

If Ai,A 2 are complex conjugate one of the eigenvalues entirely determines the other, 
and the joint probability density automatically reduces to the one-level density for 
complex eigenvalues 


^C(A) = ^IMA)l e-2(Re(A)Mm(A)^)erfc(2|im(A)|). 


(14) 


The factor of two accountes for the fact that the eigenvalues are indistinguishable. 
The total one-level density of a complex eigenvalue is then given by 

Ri(A) = i?®(Re(A),Im(A)) +5(Im(A))i?f(Re(A)). (15) 

Figure [^depicts the numerically obtained histogram of the eigenvalues of a large 
number of split-complex Hermitian matrices sampled from the distribution (§. We 
also plot the probabilities conditional on real and complex eigenvalues respectively 
in comparison with the analytical results. We observe a good agreement between 
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X 


Figure 1. One-level density for the 2x2 split-complex Hermitian ensemble 
§. The top left panel shows the histogram of the eigenvalues using numerically 
obtained eigenvalues from 200,000 matrices. The right panels show the histogram 
of the conditionally complex eigenvalues (top) in comparison to the analytical 
result (bottom). The bottom left panel shows the numerical and analytical 
distribution for conditionally real eigenvalues. 




the numerics and the analytical results; in particular, there is a sharp peak in the 
histogram for Im(A) = 0, reflecting the 5-function in the analytic distribution. 

Note that the distributions (131 and (141 are similar to those obtained from 

In fact, the two ensembles can be explicitly 


37 


the 2x2 real Ginibre ensemble 
related. The 2x2 split-complex Hermitian ensemble can be brought, via a parameter- 
independent orthogonal transformation on the 4x4 matrix representation, to a matrix 
which is spectrally equivalent to two copies of a matrix from the 2x2 real Ginibre 


ensemble, with elements distributed as A/'(0, 
that 


O' 


( 

0 

6 

7 \ 


( 

a 

b 

0 


0 

Ai 

7 

5 

0 = 


d 

c 

0 

0 

5 

-7 

A 2 

0 



0 

0 

a 

d 

V -7 

5 

0 

A 2 ) 


\ 

0 

0 

b 

c / 


\) instead of 1). Explicitly we find 


(16) 


with 




1 

0 

\ 

0 = ^ 
V2 

0 

1 

0 

1 

1 

0 

-1 

0 


11 

0 

1 

0 / 

M 

1 

= 5 + 

7 , 

and c 

= Ai. 


( 17 ) 


identically distributed as A/*(0, ^). A similar relation holds in the N x N case. 
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Let us now consider the split-quaternionic case. Let 

^ ^ (5 - i/r - jq - kcr 

i5 + i/r + jq + kcr h-2 

with Ai 2 , d, & R; that is, we have the probability distribution 

P(Ai 2, d, At, 7, cr) = 

Transforming to the variables Ai ^25 <5? cind 7 , we obtain 

32e-2('^?+'^2+4(7''+a-''))|y^ - A 2 I 


]P(Ai,2,d,At,7:0-) = 


(18) 


(19) 


( 20 ) 


^V(Ai - A2)2 + 4n2 
where = ^2 _|_^2 _ ^2 _ proceed to integrate over the regions where the 

original parameters remain real. To ensure that Ai 2 € IR we need —h > 0. 

Rearranging this inequality as a limit on At and integrating over fx, followed by 
integrating over S over the region — 1 _ ^2 _|_ u 2 ^ leaves 

P(Ai.2,7,a) = 32|Ai-A2| ^_2(,^+,2+4(,-+..))g^^^ ^21) 

TT^ 

where f (Ai, 2 )^ = ^'^^-^ 2 ) — 1 _ ^2 _|_ ^ 2 _ ensure that (5 S IR, we additionally require 
'^^(Ai, 2 )^ > 0. Now, transforming our variables ( 7 , 0 -) into polar coordinates {r,(j)) and 
integrating over (j) G [0, 27r] leaves the joint probability density function of Ai, A 2 , 


P(A 


^^A2)^6^A^7^e-2(A?+A=)7^(A,_2), 


( 22 ) 


where 72.(Ai^2) = /r dr. The region T is different for real 

eigenvalues in comparison to complex conjugate eigenvalues. 

In the case of real eigenvalues, the reality of S implies that we integrate over 
r G [0,oo]. We thus find the joint probability distribution for real eigenvalues as 

P7Ai,A2)= ^(Ai-A2)V2("?+"^) 

TT 

+ l^i^2l e-4AiA2ei.fc(72|Ai-A2|). (23) 


Integrating over one of the eigenvalues we find the one-level density in the case of real 
eigenvalues as 


<(A) = 







(24) 


In the split-quaternionic case, the probability of obtaining real eigenvalues is 1 — , 

slightly smaller than in the split-complex case. 

If Ai and A 2 are complex conjugate, however, the reality of S means that 
7 ^ -I- > Im(Ai^ 2 )^- That is, we now have to integrate over r G [|Im(Ai_ 2 )|, c»]. 

Accounting for the fact that the eigenvalues are indistinguishable, we thus obtain the 
one-level density for complex eigenvalues 


R®(A) = 2Y^|Im(A)|e"'‘(^^®(^^^'+^^“(^^^"). (25) 

Again we can express the total one-level density of a complex eigenvalue in the form 

(pl. 
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Figure 2. As figure^ for the split-quaternionic case. 


Figure [^depicts numerically obtained eigenvalue densities for a large number of 
matrices sampled from the distribution in comparison with the analytical results. 
We observe a dip in the one-level density for real eigenvalues, similar to the one in 
the GUE one level density for 2x2 matrices. The exact form of the distribution (24) 
however, is distinctively different from the GUE case. 

Let us finally derive the level spacing distributions in the case of real eigenvalues 
for the split-Hermitian ensembles. These are the distributions of the spacings 
s oc |Ai — A 2 I, subject to the requirement that the mean level spacing is one. For 
this purpose we subsitute /ii = Ai — A 2 , /i 2 = Ai -I- A 2 , in the joint probabilities of the 
real eigenvalues. For the split-complex case this yields 

1 


4v7r 


(26) 


Integrating /i 2 over the whole real line, normalising and rescaling to ensure that the 
mean spacing is one, we obtain 


(s) = ^se- 


(27) 


This is identical to the level spacing for the GOE, as has been noted before for the 
real Ginibre ensemble 


26 


For the split-quaternionic Hermitian ensemble we again substitute /ii and ^2 hr 
the joint probability distribution for real eigenvalues to obtain 

^ - l/ii 


p(m 


i,2j — —e 

TT 


Ml M2 


2 v^ 




(28) 
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Figure 3. Level spacing distributions for real eigenvalues from 100,000 matrices 
from the split-complex (left) and split-quaternionic (right) ensembles. 


After integrating over and appropriate rescaling we obtain 


s = 


/2a3 ^e- 


+ ^se“^"erfc(v^s)^ 


2 ^ 2-1 


(29) 


1 / 3V2-sinh“bl) 

“ = I i2V2-l)V^ 


VO 

Figure ([^ depicts the numerically obtained histogram of the level spacings of 
the eigenvalues of a large number of split-complex Hermitian and split-quaternionic 
Hermitian matrices sampled from distributions ([^ and (|^. We observe an excellent 
agreement with the analytical reasults. Note that the level repulsion is linear for both 
ensembles. 

In summary, we have shown that the space of PT-symmetric matrices is 
isomorphic to the space of split-Hermitian matrices. We have introduced two new 
random matrix ensembles of split-complex and split-quaternionic Hermitian Gaussian 
random matrices, and derived analytic expressions for the 2x2 case. We have 
demonstrated a relation between the real Ginibre ensemble and the split-complex 
Hermitian ensemble. It is a challenging but achievable task to derive analytic results 
for general N x N matrices for the new ensembles. The split-complex Hermitian 
ensemble might help in understanding some of the still-unknown spectral properties 
of the N X N real Ginibre ensemble. We have conjectured that the newly derived 
ensembles constitute universality classes for PT-symmetric systems. Hitherto only few 
PT-symmetric model systems have been investigated with respect to their statistical 
spectral properties. In 38 , GOE-type statistics were observed for a PP-symmetric 




Dicke model in the special case of real eigenvalues, which is consistent with the spectral 
properties of the split-complex ensemble introduced here. It is an important task to 
compare the spectral features of the new ensembles to those of further PP-symmetric 
model systems. 
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